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Abstract 

The Lagrange-mesh method has the simplicity of a calculation on a mesh 
and can have the accuracy of a variational method. It is applied to the study 
of a confined helium atom. Two types of confinement are considered. Soft con¬ 
finements by potentials are studied in perimetric coordinates. Hard confinement 
in impenetrable spherical cavities is studied in a system of rescaled perimetric 
coordinates varying in [0,1] intervals. Energies and mean values of the distances 
between electrons and between an electron and the helium nucleus are calculated. 
A high accuracy of 11 to 15 significant figures is obtained with small computing 
times. Pressures acting on the confined atom are also computed. For sphere radii 
smaller than 1, their relative accuracies are better than 10“^^. For larger radii 
up to 10, they progressively decrease to 10“^, still improving the best literature 
results. 


1 Introduction 

Many numerical techniques exist for quantum-mechanical calculations in conhguration 
space. Among them, two main qualities may be searched: accuracy and simplicity. 
However, they are not often encountered simultaneously. For some problems, the 
Lagrange-mesh method has the accuracy of a variational method and the simplicity of 
a calculation on a mesh [B m m SI i5]. This approximate variational method involves 
a basis of Lagrange functions, i.e. inhnitely differentiable functions vanishing at all 
mesh points of a Gauss quadrature, except one. With the help of the associated Gauss 
quadrature, all matrix elements are very simple. In particular, the matrix elements 
of the potential are approximated by values of the potential at mesh points like in 
collocation methods. 

The striking property of the Lagrange-mesh method is that, in spite of its simplicity, 
it has essentially the same accuracy as a variational calculation performed with the 
same Lagrange basis. This performance is not well understood yet [3]. However, the 
accuracy of the method depends on the validity of the Gauss quadrature. Hence, the 
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method can be very bad in the presence of singnlarities of the potential. This problem 
can sometimes be cured by a so-called regularization ElEliaE]. The method has been 
successfully applied to many problems in atomic, molecular and nuclear physics (see 
Ref. [S] for a review). It is particularly useful to solve coupled-channel problems in the 
continuum EE] or three-body problems mu- 

The aim of the present paper is to apply the Lagrange-mesh method to a three- 
body problem: a helium atom conhned in some environment [Ini[lll[l2l[l3l[ll[l5l[l6]. 
The conhnement can simulate a helium gas under pressure or helium atoms trapped in 
some molecule, carbon cluster or crystal. As we show below, existing techniques na 
are very convenient for soft conhnements by potentials but can not be applied for an 
atom conhned in an impenetrable spherical cavity. The difficulty arises from simul¬ 
taneously meeting two different constraints: regularizing the singularities of the three 
Coulomb terms and forcing the conhnement. The former condition is easily treated 
in perimetric coordinates usiEn] which automatically regularize the singularities but 
which necessarily extend over the whole conhguration space. Conhnement is easily 
and accurately treated on a Lagrange mesh over a hnite interval for hydrogen |20l|5]. 
Here we show that a new coordinate system can be built for helium which keeps the 
regularization of the potential but over an impenetrable spherical cavity. 

The principle of the Lagrange-mesh method is recalled in section HJ In section |31 
the problem of the conhned helium atom is presented. In section 01 the Lagrange- 
mesh method in perimetric coordinates is recalled. A new system of coordinates is 
introduced in section |5] as well as its Lagrange-mesh implementation. Results for the 
diherent types of conhnement are discussed on section O Concluding remarks are 
presented in section [71 

2 Principle of the Lagrange-mesh method 

Let us consider N mesh points Xi associated with a Gauss-quadrature approximation 




( 1 ) 


k=l 


The weight coefficients Aj are also called Christohel numbers. Lagrange functions are a 
set of N orthonormal functions fj {x) associated with this mesh verifying two conditions 
ID HIS]- (i) They satisfy the Lagrange property 



( 2 ) 


i.e., they vanish at all mesh points, but one. (ii) The Gauss quadrature is exact for 
products of two Lagrange functions. 


Let us consider a particle of mass m in a potential V (x). This basis is used in an 
approximate variational calculation with the trial function 


N 


^(x) =Y,Cjfj{x). 


(3) 
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The matrix elements of potential V(x) are calculated at the Gauss approximation as 

N 

{fi\y\fj) ~ ^kfi{xk)V{xk)fj{xk) = V{xi)5ij (4) 

k=l 

because of the Lagrange property ([2]). The variational equations then take the form of 
mesh equations BSIE] 

^ f \ 

Y. G = Eci, (5) 

where the exact or approximate matrix elements T^- = (/j| — (P/dx‘^\fj) have simple 
known expressions as a function of the zeros Xj and Xj. See Refs. BIlllS] for details. 

This simple approximation can be very accurate with small numbers of mesh points 
when the Gauss approximation is valid for the potential matrix elements, i.e. when the 
potential and its derivatives have no singularities. In the presence of singularities, 
the basis can sometimes be regularized by multiplying the Lagrange functions by a 
convenient factor R{x) [2llll|5], 

The regularized functions fj{x) still verify the Lagrange conditions ([2]) but they are 
not orthogonal anymore, in general. However, they are still orthogonal at the Gauss- 
quadrature approximation and can still be treated as orthonormal in the Lagrange-mesh 
method without signihcant loss of accuracy mn- Such a regularization is also useful 
when the particle is conhned by an impenetrable wall |20] (see section [5]). 

3 Confined helium 

We consider a two-electron atom with an inhnite-mass nucleus of charge Ze. This 
nucleus is hxed and the electrons are characterized by coordinates ri and r 2 with 
respect to this nucleus. In atomic units h = me = oq = e = 1, where me is the electron 
mass and Oq is the Bohr radius, the Hamiltonian of the helium atom reads 

// = r + V'c = -iAi-lA2-(7) 

2 2 ri r2 ri2 

where 

ri2 = ri-r2 (8) 

and Ai and A 2 are the Laplacians with respect to Vi and r 2 . 

For a free atom, the wave functions of the bound states must vanish at inhnity. A 
conhnement can be introduced in the problem either by forcing the wave function into 
some spherical cavity (hard conhnement) or by adding a conhning potential to H (soft 
conhnement). 

Hard conhnement is obtained with the conditions 

ri < R, r2 < R. (9) 
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The wave function ’" 2 , ^ 12 ) of an S state must thus verify 

ijiR, r2, ru) = '0(ri, R, r^) = 0. (10) 

Soft conhnement can be obtained by adding a potential Konf(^"i, ^ 2 ) which tends 
to a large positive constant or to inhnity when ri or r 2 tends to inhnity. The role of 
this potential is to reduce the probability density of presence of the electrons at large 
distances. 

Let us start with the soft conhnement since it can be treated with the same code 
as for the free atom with only a tiny modihcation m- 

4 Lagrange mesh for soft confinement 

4.1 Perimetric coordinates 

The system of perimetric coordinates [IHl US] is very convenient for Lagrange-mesh 
calculations of three-body systems because the three dimensioned coordinates are in¬ 
dependent from each other and vary from zero to inhnity. Moreover, the volume element 
automatically regularizes the singularities of the three Coulomb potentials HU. 

The perimetric coordinates are composed of three Euler angles and the three coor¬ 
dinates 


X = n - r2 + ri2, 

y = -n + r 2 + ri 2 , 

z = ri + r2-ri2. ( 11 ) 

The volume element of the dimensioned coordinates reads 


dV = {x + y){y + z){z -1- x)dxdydz. 


( 12 ) 


In perimetric coordinates, the Coulomb potentials become 

2Z 2Z 2 


Vc{x,y,z) = 


z + X y + z X + y 


(13) 


With the volume element dV, the integrand in matrix elements of this potential is 
bounded everywhere. Hence the Gauss quadrature and the Lagrange-mesh method are 
accurate. 

The kinetic-energy operator T is rather complicated [22] • It is convenient to write 
its matrix elements in a symmetric form [ssiini, 

/•oo roo fOO 3 P)T? 

{F\T\G) = 2 dx dy dz (14) 

Jo Jo Jo OXi OXj 


where {xi,X 2 ,X 3 ) = {x,y,z). The coefficients Aij are given by 


All = x{y + z){x + y + z) xz{z + x), 

^22 = yz{y + z) + y{z + x){x + y + z), 

^33 = yz{y + z) + xz{z + x), 

Ai2 = ^421 = 0 , 

Ai 3 = A31 = -xz{z + x), 

^23 = ^32 = —yz{y + z). ( 15 ) 
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4.2 Lagrange mesh and functions 

The Lagrange-Laguerre functions are defined as [1] 

fj(^) = (-lyxf e-"'^ (16) 

where L]s[{x) is the Laguerre polynomial of degree N and the Xi are its zeros, 


Ln(Xi) = 0 . 


(17) 


Notice that the Lagrange functions are linearly independent polynomials of degree 
— 1 multiplied by an exponential which is the square root of the Laguerre weight 
function exp(—x). The basis is thus equivalent, for example, to a basis formed of the 
Laguerre polynomials of degrees 0 to — 1 multiplied by exp(—x/2). 

The functions fj{x) are associated with the Gauss-Laguerre quadrature IS]. They 
verify the Lagrange conditions ([2]) and integrals of products of two Lagrange functions 
are exactly given by the Gauss quadrature since the integrand is the product of the 
Laguerre weight function exp(—x) by a polynomial of degree 2N — 2 [2l]. Functions 
f[TB]l are exactly orthonormal over (0, oo). 

The hrst derivative of a one-dimensional Lagrange-Laguerre function at mesh points 
is given by 


for i ^ j and by 






1 

2xj' 


(18) 


(19) 


Three-dimensional mesh and basis are obtained as follows. Let Xp {p = 1,, N^), 
Pq (g = 1,..., Ny) and Zr {r = ..., N^) be the zeros of Laguerre polynomials with 

respective degrees Ny and N^. Three-dimensional Lagrange functions Fijk{x,y, z) 
associated with the mesh {h^Xp, hyPy, hzZr) are dehned by 


F,jk(x, y, z) = Ni.l''‘fP’\xlh^)ff'\ylhy)f'P‘\zlh, 


( 20 ) 


The functions are given by expression flTB]) with N replaced by Ny or N^. The 
corresponding Ghristoffel numbers are denoted as Aj, pj and z/fc. Scale parameters 
hy and are introduced in order to £t the different meshes to the size of the actual 
physical problem. The normalization factor Mijk is dehned as 


•^ijk hxhyh^i^hxXi “1“ hyPj^i^hxXi h^Zk)(yhyyj -I- h^z^). (^^) 

The Lagrange functions Fijk{x,y, z) satisfy the Lagrange property 

Fijk{hx^pt hyVqi hzZj.) (-A/ijfcAj/XjX'^) ^ ^jq^kr 5 (22) 


i.e., they vanish at all points of the three-dimensional mesh, but one. With the volume 
element o, they are not orthogonal but the scalar product is calculated 
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with the Gauss-quadrature approximation as duidjjidkk'- They are thus treated as an 
orthonormal basis in the method. 

The kinetic-energy matrix elements are given by 


-2 

X 


{Fiijik'\T\Fijk) ~ Khyh^ 

Aii{h:,Xn, hyl/j, hzZk)XnfY"^'{Xn)fY"^'{Xn)h., 

n 

F^iiAkk’ Y A22iKXi, hyyn, h^Zk)ynfY'‘^\yn)f\^''^\yn)hy‘^ 

n 

F^uAjj' Y ^SsiKXi, hyVj, h^Zn)iynfY^^'{Zn)flY'{Zn)h~‘^ 

n 

+Skk' [AuiKxi, hyi/j/, h^Zk){XiiJ,jY‘^ {yj') 

+Ai2{Kxi>, hyyj, h^Zk){Xi>iJ,jY^‘^fY"'''{xii)f^!^^'''{yj) {Khy)~^ 
+5jf Ai3{h^Xi, hyyj, h^Zk'){XiUk>Y^‘^fY"'''{xi)fY"^'{zk>) 
+Ai3{Kxi>, hyyj, h,Zk){XeiZkY^‘^fY"'^'{xiYfiY\zk)] {KhY~^ 




A23{Kxi, hyyj, h^Zk'){yjiZk'y^^ff'''’'{yj)fY^^'{zk') 


+A23iKxi,hyyj>,h^Zk)i/J-fiykY^‘^fj^'''''iyj')fiY\zk) {kyh^ (23) 


From now on, we consider = Ny = N which implies Xi = yi, Xi = yi, and = 

hy = h. 

The Lagrange functions are used as a variational basis to expand an S'-wave trial 
function, 


N i-a 

'4^{.x,y,z) = YYY Cijk[2{l + 5p)] ^/‘^[Fijk^x, y, z) ± Fjik{x, y, z)], (24) 

i=l j=l k=l 

where a = 0 in the symmetric case and 1 in the antisymmetric case, and j < i — (t 
because of the symmetry with respect to the exchange of electrons 1 and 2. The matrix 
representing the total potential V = Vc F Konf is immediately obtained with a triple 
Gauss quadrature and the Lagrange property fl22l) as 


i^F^tjik^\\z\Fj^jiY ~ F(^hx{^hxj^h^Zk')SiifSjj/Skk'‘ (^^) 

The variational calculation then reduces to mesh-like equations 

N i—a Nz 

Ei; i:{(l + hr'F + hi')-F(F.'j't'\T\F.jt) ± {Fff,.\T\F,it)] 

i=lj=\k=l 

“[“[V(/iX2, hXj^ ^zZk) 

Energies and wave functions are obtained from the eigenvalues and eigenvectors of a 
large sparse symmetric matrix. 

Mean values of a multiplicative operator 0{x,y,z) are simply given at the Gauss 
approximation by 

N i-a Nz 

{'ip\0{x,y,z)\Y) ^ I YYY Cijk[Oihxi, hxj, hzZk) + 0{hxj, hxi, h^Zk)]. (27) 

i=lj=lk=l 
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5 Lagrange mesh for hard confinement 


5.1 Rescaled perimetric coordinates 

In an impenetrable spherical cavity of radius R, the perimetric coordinates are con¬ 
strained by 


0 < X < 2R — z, 

0 < 1/ < 2R — z, 

Q<z<2R. (28) 


Hence the kinetic matrix elements o must be rewritten as 


|■2R |•2R-z r2R-z 

{F\T\G) = 2l dzl dxl 9 , 

-1 


(29) 


The z dependence of the upper bounds of the x and y coordinates is a difficulty for 
the use of the Lagrange-mesh method. 

We thus introduce a new set of coordinates (u, v, w) = {ui,U 2 , Us) defined over [0,1] 
by 


Their upper values are indeed 


X 


u = 


V = 


2R-z' 

y 

2R-z' 

z 


w = 


2R' 


u = l (ri = R), 
v = l (ra = R), 
w = 1 {ri = r 2 = R, ri 2 = 0). 


Inversely, one has 


X = 2Ru{l — w), 
y = 2Rv{l — ta), 
z = 2Rw. 


(30) 


(31) 


(32) 


The volume element then becomes 


dV = {2RY{u -\- v){u -\- w — uw){v -\- w — vw){l — w)^dudvdw. (33) 

The Coulomb potential reads 

Vc{u,v,w) = ^ 


+ 


u + w — UW V + w — vw {u + v){l — w) 


(34) 
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The integrands in its matrix elements are automatically regularized by the volume 
element. The kinetic matrix elements become 


where 


with 



dF dG 

dui duj 


Bii = u{l — w) ttyb + (1 — uYwau 
B 22 = '1^(1 — w) (1 — vYwtty + aj) 

S33 = w{l - wf [vtty + UQu] , 

Bi 2 = B21 = uvw{l - w) [(r> - l)a^ + {u - l)aj , 
Bis = B31 = uw{l - wf [vay + {u- l)aj , 

B23 = Bs 2 = vw{l - wf [(n - l)a^ + Ma„] 


ttu = u + W — UW, 
tty = V + W — VW, 

b = u + v + w — uw — VW + uvw. 


(35) 


(36) 


(37) 


5.2 Lagrange mesh and functions 

Let us introduce a convenient basis over the [0,1] interval. Regularized Lagrange- 
Legendre functions are defined by 


fj{u) = (-1) 


N-j 


I Uj Pn{2u — 1 ) 

1 — Uj U — Uj 


(l-tl), 


(38) 


where the mesh points Ui are the zeros of the shifted Legendre polynomial of degree 

N, 


PN{2ui - 1) = 0. (39) 

They are associated with the Gauss-Legendre quadrature on the [0,1] interval. The 
regularized Lagrange functions (l38|) correspond to standard Lagrange-Legendre func¬ 
tions [T], shifted [25] and multiplied by the factor (1 — m)/(1 — Uj), so that they vanish 
at M = 1. In Ref. [25] on the contrary, they are multiplied by u/uj so that they vanish 
at M = 0. In Ref. ra. they are multiplied by both factors simultaneously and vanish at 
M = 0 and 1. Notice that all these types of regularized functions satisfy the Lagrange 
conditions (|2]) but are not orthogonal. Nevertheless they can be treated as orthonormal 
in Lagrange-mesh calculations without signihcant loss of accuracy 

The hrst derivative of a Lagrange-Legendre function fl38|l at mesh points is given 
by 


Uj{l 



Uj) 1 

Uj) Ui — Uj 


= (-1)'*' 


(40) 












for i ^ j and by 




2'Uj(l Mj) 


( 41 ) 


Over a three-dimensional mesh {up, Vg, Wr), where Up, Vg, are solntions of Eq. 
with possibly different valnes Nu, Ny, N^j of N, three-dimensional Lagrange fnnctions 
Fijk{u,v,w) are dehned by 




with 


Afijk = {2Rf{ui + Vj){ui + Wk- UiWk){vj +Wk- VjWk){l - Wkf. 
These functions satisfy the Lagrange conditions 

Fijki^Up^ Vq^ Wr) ^ ^ip^jq^kvi 


(42) 


(43) 


(44) 


where A*, /ij and Vk are now the weights of the Gauss-Legendre quadrature over the 
[0,1] interval with Ny and points, respectively. They are not orthogonal but 
they are orthonormal at the Gauss-quadrature approximation and are treated as an 
orthonormal basis in the method, 

i^Fi'jik'\Fijk') y Sii'6jj'5kk' ■ 

The kinetic matrix elements are given by 

« 2(2F)W.7,:^Wy;''" 

X [Sy'^kk' ^ Fii('((„, Vy Wk)KF'‘^\nn)fl^'''\Un) 

n 

F^ii'hk' y B22{Ui, Vn, Wk)finfj^"^'{Vn)fj!^"^'{Vn) 
n 

y ^ 33 (^ 4 , 

n 

+hk' [Suiui, Vf, 

+Sjj> [Bi^iui, Vj, Wk'){Xii^k'Y^‘^f!:'^"^\ui)fjf'^'’'{wk') 

+ Bi3{Ui>,Vj,Wk){Xi'l^ky^‘^fi^"^'{Ui>)flY^"’^'{Wk) 

+5ii' B23{Ui, Vj, Wk'){fijl^k'Y^'^fjY''''\vj)fi^"^\wk') 

From now on, we take Ny = Ny = N. 

An S'-wave trial function is expanded as 

N i — a Nw 

'Ip{u,v,w) = y y y ^^^*,[ 2(1 Sij)]~^^‘^[Fijk{u,v,w) ± Fjik{u,v,w)], (47) 

2=1 j=l k=l 
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where a is defined as above. The potential matrix elements simply read 




1 

X — 

R 


+ 


Ui T UiWf^ Vjj T ta/;; UjWj^ T ra/j;) 

The Lagrange-mesh eqnations are 

N i — a Nw 

E E E{(1 + ± {F,,,,,\T\F,,k)] 

i=lj=lk=l 


( 48 ) 


(49) 


Notice that obtaining energies does not reqnire calculating eigenvalues of a generalized 
eigenvalue problem since the basis is treated as orthonormal. 

Mean values of a multiplicative operator 0{u,v,w) are given at the Gauss approx¬ 
imation by 

N i—u Nw 

{'ip\0{u,v,w)\'ip) ^ i EEE Flk[0{ui, Uj, Wk) + 0{uj, Ui, Wk)]. (50) 

i=lj=lk=l 


6 Results 

All results are presented in atomic units except in Table [71 An error of a few units 
may affect the last displayed digit. 


6.1 Soft confinement 

The study of soft conhnement is performed in perimetric coordinates using the code 
of Ref. [IT] with two modihcations. (i) A conhning potential Konf is added to the 
Coulomb potential; this is a very simple addition in the code since no matrix elements 
are needed, (ii) The search of the lowest eigenvalues of the rather scarce Hamiltonian 
matrix is now performed with the Jacobi-Davidson algorithm [26]. Mean values are 
calculated with Eq. fl27|) for the distance ri 2 between electrons and for the distances ri 
and r 2 between the electrons and the nucleus. 

The calculations are performed to provide 12-15 signihcant hgures in the fastest way, 
when such an accuracy can be reached with reasonable basis sizes. The calculation 
starts with a search for an optimal domain of the scale parameters h and h^- The 
signihcant digits are obtained by comparison between several calculations with different 
numbers of mesh points. The simplest of the calculations giving the requested accuracy 
is then kept. 

As a hrst penetrable conhnement, we choose like in Ref. [12] the harmonic potential 

Rconf(ri, r2) = + r^), (51) 

where a; is a parameter. The convergence for a; = 1 is studied in Table [Has a function 
of the numbers N and N^. The total mesh size Nt = \N[N + 1)A'^ is also given. 

The parameters h and hz are roughly optimized. One observes a fast convergence 
for h = 0.15 — 0.20 and hz = 0.15 — 0.20 where very accurate results are already obtained 
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Table 1: Convergence of the gronnd-state energy and the mean interparticle distances 
of a helium atom conhned by the harmonic potential fl5T]) with a; = 1 as a function of 
the numbers N and Nz of mesh points and the total mesh size Nt- 


N 

Nz 

Nt 

h 

hz 

E 

{ri2) 

(ri) = (r 2 ) 

10 

10 

550 

0.15 

0.20 

-2.073 075 9 

1.085 8179 

0.723 715 2 




0.20 

0.20 

-2.073 035 3876 

1.085 686 26707 

0.723 644 4069 

15 

15 

1800 

0.15 

0.20 

-2.073 035 362 032 5 

1.085 685 768 80716 

0.723 644141796 32 




0.20 

0.20 

-2.073 035 362 0477 

1.085 685 768 60892 

0.723 644141692 94 

20 

20 

4200 

0.15 

0.20 

-2.073 035 362 05189 

1.085 685 768 624 21 

0.723 644141700 93 




0.20 

0.20 

-2.073 035 362 05154 

1.085 685 768 624 56 

0.723 644 14170114 

25 

20 

6500 

0.15 

0.20 

-2.073 035 362 05195 

1.085 685 768 624 25 

0.723 644 141700 96 




0.20 

0.20 

-2.073 035 362 05188 

1.085 685 768 624 25 

0.723 644 141700 96 

25 

25 

8125 

0.15 

0.20 

-2.073 035 362 05194 

1.085 685 768 62419 

0.723 644 141700 92 




0.20 

0.20 

-2.073 035 362 05186 

1.085 685 768 62399 

0.723 644141700 80 

30 

25 

11625 

0.15 

0.20 

-2.073 035 362 05194 

1.085 685 768 62413 

0.723 644141700 88 




0.20 

0.20 

-2.073 035 362 05194 

1.085 685 768 624 22 

0.723 644141700 93 

30 

30 

13950 

0.15 

0.20 

-2.073 035 362 05193 

1.085 685 768 624 21 

0.723 644141700 91 




0.20 

0.20 

-2.073 035 362 05191 

1.085 685 768 624 23 

0.723 644 141700 93 

Ref. 

m 




-2.073 035 




Table 2: Ground-state energy and mean interparticle distances of a helium atom con¬ 
hned by the harmonic potential flSTj) as a function of the numbers N and Nz of mesh 
points and the scale parameters h and 


UJ 

N 

Nz 

h 

hz 

E 

(^ 12 ) 

(^ 1 ) = {^2) 

0 

30 

25 

0.30 

0.35 

-2.903 724 377 034 0 

1.422 070 255 565 9 

0.929472 294 873 7 

0.01 

30 

25 

0.30 

0.35 

-2.903 605 041422 9 

1.421940 016 512 7 

0.929 395 600 292 0 

0.05 

30 

25 

0.30 

0.35 

-2.900 748 507134 5 

1.418 859 663 654 5 

0.927580 543426 0 

0.1 

25 

25 

0.25 

0.30 

-2.891910 703103 4 

1.409 741024 0940 

0.922195 532 048 2 

0.25 

25 

25 

0.25 

0.25 

-2.833 069 315 3970 

1.359 681803 628 3 

0.892 361239 0213 

0.5 

25 

20 

0.20 

0.20 

-2.648 703149 419 7 

1.256 367 735 297 2 

0.829 664522 705 9 

1 

25 

20 

0.20 

0.20 

-2.073 035 362 0519 

1.085 685 768 624 2 

0.723 6441417010 

2 

20 

20 

0.15 

0.15 

-0.493173 861504 2 

0.879 296 076448 4 

0.592 319 993 998 0 

5 

20 

20 

0.10 

0.10 

5.555 021418 896 5 

0.6193775846175 

0.422 804512 778 6 

10 

20 

20 

0.06 

0.06 

17.1621913740574 

0.459 421428 2815 

0.316 272 749 895 5 


for N = Nz = 20. The results are quite insensitive to variations of hz in this interval. 
Rounding errors start to slightly deteriorate the accuracy around N = Nz = 30. For 
small values of u, the convergence is very similar to the free-atom case. See e.g. Table 1 
of Ref. [1^ . 

The Lagrange-mesh results for various u values are presented in Table [2l The values 
for the free atom {uj = 0) have an accuracy of about 10“^^ like in Ref. [1^. A relative 
accuracy of 10“^^ is already obtained with N = Nz = 20. The conhnement reduces 
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Table 3: Convergence of the gronnd-state energy and the mean interparticle distances 
of a helium atom conhned by potential with Vq = 100 and R = 4. The scaling 
parameters are h = 0.10 and hz = 0.15. 


N 

Nz 

Nt 

E 

(n2) 

(ri) = (ra) 

10 

10 

950 

2.3895176 

0.717239 6 

0.487082 6 

15 

15 

1800 

2.389 521083 273 

0.717240 933 789 

0.487083 296 959 

20 

20 

4200 

2.389 521083 37013 

0.717240 933 82109 

0.487083 296 966 05 

25 

25 

8125 

2.389 521083 37011 

0.717240 933 82111 

0.487083 296 966 08 

Ref. [12] 


2.389 531 




the size of the atom when u: increases while increasing its energy. The scale factors for 
the weakest conhnement are inspired by those of a free atom. When u increases, they 
progressively decrease. To obtain a constant accuracy (thirteen signihcant digits), the 
numbers N and Nz can both decrease when the conhnement becomes stonger. The 
present energies conhrm the six-digit results of Ref. [12] for w < 1, up to a possible 
rounding. 

Table 4: Ground-state energy and mean interparticle distances of a helium atom con¬ 
hned by potential (l52|) as a function of the numbers N and Nz of mesh points and the 
scale parameters h and hz- 


R 

N 

Nz 

h 

hz 

E 

(^ 12 ) 

{ri) = {^ 2 ) 

Vo = 

: 25 







1 

20 

20 

0.10 

0.10 

9.005 664919 615 0 

0.574015409 803 2 

0.392137172 502 8 

2 

20 

20 

0.10 

0.10 

2.205 633 821802 6 

0.732 574434299 6 

0.496 887442 6881 

3 

20 

20 

0.10 

0.15 

0.058 610162109 7 

0.843146 864040 8 

0.568 830 018 0991 

4 

20 

20 

0.10 

0.15 

-0.937974 874 2074 

0.926139 670492 3 

0.622 265 713 845 8 

5 

20 

20 

0.15 

0.15 

-1.493 463 026 587 6 

0.991265 922 0321 

0.663 853 999 5004 

10 

20 

20 

0.20 

0.30 

-2.440 318186 433 3 

1.180598 9951242 

0.782 896 512 886 3 

25 

25 

20 

0.30 

0.35 

-2.814 574 813 589 3 

1.346 698 806 392 7 

0.884 554 240409 3 

100 

25 

20 

0.30 

0.35 

-2.897 788 906 088 7 

1.415 741535 9464 

0.925 741063 258 9 

Ro = 

: 100 







1 

25 

25 

0.08 

0.08 

25.983 328 462164 7 

0.407916 008 582 9 

0.281363150 279 2 

2 

20 

20 

0.10 

0.10 

9.919 341323 7201 

0.543 627679 286 8 

0.372 439 8341115 

3 

20 

20 

0.10 

0.10 

4.812 589 834 5331 

0.640 875449 3044 

0.436 905 248 6715 

4 

20 

20 

0.10 

0.15 

2.389 521083 3701 

0.717240 933 8211 

0.487083 296 9661 

5 

20 

20 

0.10 

0.15 

1.006 687 798 027 9 

0.780 044 762 396 9 

0.528 060 835 913 2 

10 

25 

20 

0.20 

0.25 

-1.474872 423 516 3 

0.985 908 228 888 3 

0.660 518 629 618 2 

25 

25 

20 

0.30 

0.35 

-2.586 842 706183 4 

1.230 8212121751 

0.813 965 678 584 9 

100 

25 

20 

0.30 

0.35 

-2.880 323 145 5112 

1.398 590 911428 2 

0.915 588 071052 0 


As a second potential with a penetrable conhnement, we choose like in Ref. [12] 

Rconf(ri, ra) = W (2 - , (52) 
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where Vq and R are parameters. 

The convergence is stndied in Table [3] for Vq = 100 and i? = 4 as a function of N 
and Nz- Good values of the scale factors are much smaller than for the free atom na. 
Here they are roughly optimized as h = 0.1, hz = 0.15. The convergence is very fast. 
Good results are already obtained with N = Nz = 10. An accuracy of about 14 digits 
is obtained with N = Nz = 20 for the energy and for the mean distances. 

The Lagrange-mesh energies and mean distances are displayed in TablelUfor Vq = 25 
and 100 and some values of R. The scale parameters are rather small for R = 1 
and must increase with R. The ground-state energy is positive for the smaller R 
values but the three-body system is nevertheless deeply bound with respect to the 
asymptotic value 2Vo of the conhnement potential. Hence the wave functions decrease 
rather fast. Even for R = 10, the properties signihcantly differ from the free atom. 
The mean distances are scaled down with respect to the free atom im but (ri 2 ) and 
(ri) indicate that the conhnement does not change much the shape of the electron 
distribution. Indeed, the ratio (ri 2 )/(ri) progressively increases from the value 1.45 
at strong conhnement (i? = 1, Vq = 100) to the free value 1.53. Our results conhrm 
the four- and six-digit energies of Ref. [12] except for their last two digits displayed at 
R = A and 5. 

6.2 Hard confinement 

Now we consider a nucleus conhned in an impenetrable sphere. There are no scale 
parameters as the size is hxed by the radius R. Otherwise, the code for the new 
coordinate system follows the same philosophy as the previous one. The total mesh 
size is Nt = ^N{N + l)N^i; for singlet states and Nt = ^N{N — l)Nw for triplet states. 
The mean values are calculated from (l50D with (l3^ . 

The convergence of Lagrange-mesh results for the ground state is studied in Table [5] 
for various values of R, as a function of N and N.^. For large R, the convergence is 
rather slow but it is extremely fast for small R. One obtains 15 signihcant hgures 
for R = 0.1 and 1 but only 11 for R = 10. The situation is reversed with respect to 
soft-conhnement calculations in perimetric coordinates where larger bases are needed 
at small R. The results are compared with Refs. [ini HU El ng. A good agreement 
with six signihcant digits is observed for the most recent results. 

Energies and mean distances for the ground state are gathered in Table [HI for a 
number of R values. The energies conhrm the six-digit variational energies of Ref. [T2] 
from 0.5 to 10 and Ref. [16] from 1 to 10. The present coordinate system allows a 
better description of smaller R values. The earlier variational energies of Ref. [10] have 
an absolute accuracy below 0.001. This reference is the only one to present (^ 12 ) mean 
distances. Their accuracy increases from 10“^ to 0.002 when R varies from 0.5 to 6. 
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Table 5: Convergence of the ground-state energy and the mean interparticle distances 
of a helium atom conhned in a sphere of radius i? as a function of the numbers N and 
Njju of mesh points and the total mesh size Nt- 


N 


Nt 

E 

(^ 12 ) 

(n) = (rs) 

R = 

0.1 





8 

8 

288 

906.562 4079 

0.069 580 385 8 

0.049 5012410 

10 

10 

550 

906.562 422 907 

0.069 580 382 888 07 

0.049 501246 332 6 

10 

15 

825 

906.562 422 919 887 

0.069 580 382 88417 

0.049 501246 340 126 

15 

15 

1800 

906.562 429 919 886 

0.069 580 382 88417 

0.049 501246 340 120 

15 

20 

2400 

906.562 422 919 888 

0.069 580 382 88416 

0.049 501246 340 121 

Ref. [H] 


906.562 423 



R = 

1 





10 

10 

550 

1.015 754 975 53 

0.643 664 253 93 

0.441 796 632 210 

10 

15 

825 

1.015 754 975 90 

0.643 664 253 8710 

0.441 796 632 098 55 

15 

15 

1800 

1.015 754 976 048 33 

0.643 664 253 878 0 

0.441 796 632103 28 

15 

20 

2400 

1.015 754 976 048 41 

0.643 664 253 8779 

0.441 796 632103 31 

20 

20 

4200 

1.015 754 976 048 64 

0.643 664 253 878 0 

0.441 796 632103 34 

20 

25 

5250 

1.015 754 976 048 66 

0.643 664 253 8778 

0.441 796 632103 41 

25 

30 

9750 

1.015 754 976 048 67 

0.643 664 253 878 8 

0.441 796 632103 34 

30 

35 

16275 

1.015 754 976 048 54 

0.643 664 253 876 9 

0.441 796 632103 24 

Ref. [in! 


1.015 870 

0.643 938 


Refs. [laill] 

1.015 755 



R = 

10 





15 

15 

1800 

-2.904 79 

1.420 767962 81 

0.928 645 755 91 

15 

20 

2400 

-2.903 7272 

1.422 072 790 04 

0.929 474322 57 

20 

25 

5250 

-2.903 724 382 6 

1.422 070 196 74 

0.929 472 26315 

25 

30 

9750 

-2.903 724 375 625 

1.422 070173 66 

0.929 472 25158 

30 

30 

13950 

-2.903 724 375 691 

1.422 070172 98 

0.929 472 25122 

30 

35 

16275 

-2.903 724 375 668 

1.422 070172 34 

0.929 472 250 87 

35 

35 

22050 

-2.903 724 375 687 

1.422 070172 96 

0.929 472 25121 

Refs, [mi US] 

-2.903 724 




14 







Table 6: Ground-state energy and mean interparticle distances of a helium atom con¬ 
fined in a sphere of radius R as & function of the numbers N and N^j of mesh points. 


R 

N 


E 

(G2) 

(ri) = (r 2 ) 

0.1 

15 

15 

906.562 422 919 888 

0.069 580 382 884 2 

0.049 501246 3401 

0.2 

15 

20 

206.151712 932 762 

0.138 365 789 025 0 

0.097969 984 584 9 

0.3 

15 

20 

82.334517028118 

0.206 218 777530 

0.145 352111929 9 

0.4 

15 

20 

40.980 280 331705 1 

0.273 001759106 

0.1915916590472 

0.5 

15 

20 

22.741302 819 133 5 

0.338 577477653 

0.236 631213 255 9 

0.6 

15 

20 

13.318127241828 1 

0.402 809 660 984 

0.280412 455464 5 

0.7 

15 

20 

7.925 216 046 579 6 

0.465 563 856457 

0.322 876 819 396 5 

0.8 

15 

20 

4.610 407 554 239 4 

0.526 708460 703 

0.363 966 2772674 

0.9 

15 

20 

2.463 235 988 6210 

0.586 115 943465 

0.403 624 249 990 8 

1.0 

15 

20 

1.015 754 976 048 4 

0.643 664 253 878 

0.441796 632103 3 

2.0 

20 

20 

-2.604 038 275 176 2 

1.097202 490172 

0.733 956 380 589 2 

3.0 

25 

25 

-2.872 494 886 475 7 

1.322 925 949 994 

0.871985 220 059 

4.0 

25 

25 

-2.900 485 763 363 2 

1.399 878118648 

0.916 962 790 629 

5.0 

30 

30 

-2.903 410 849 203 1 

1.418 255199 620 

0.927363 648 349 

6.0 

30 

30 

-2.903 695 908 930 7 

1.421528 373 775 

0.929176 980 774 

7.0 

35 

35 

-2.903 7219115318 

1.422 002 693115 

0.929435 862 703 

8.0 

35 

35 

-2.903 724170 713 6 

1.422 062 567876 

0.929468 183475 

9.0 

35 

35 

-2.903 724 360 207 

1.422 069437714 

0.929471860 371 

10.0 

35 

35 

-2.903 724 375 687 

1.422 070172 936 

0.929472 251212 


From the radius dependence of the energies E{R), one can deduce the pressure 
acting on the conhned atom [ini[l2l[l6], 


1 dE 


(53) 


The derivative is performed numerically with a 4-point hnite-difference formula. The 
results are presented in Table [71 As before, the signihcant digits are estimated by 
comparison between several calculations differing by the number of mesh points. Only 
those digits are displayed in Table [71 The results are presented both in atomic units 
and in atmospheres (1 a.u. = 2.903 628 236 775 x 10® atm). 

A high accuracy is reached up to i? = 2. Beyond that value, the differences between 
neighboring energies become tiny and the relative error increases. The results are 
compared with several earlier variational calculations. The pressures of Ref. HDl have 
a relative accuracy better than 10“^ at R = 0.5 and still better than 10“® up to i? = 2. 
They become quite poor beyond i? = 4 as could be expected from the comparison 
between their variational calculations with different basis sizes. The results of Ref. [12] 
are consistent with ours for the three displayed digits except at i? = 5 and 8. The 
pressures of Ref. [IB] have a relative accuracy varying between 0.001 and 0.006. 

Energies and mean distances for the hrst excited singlet level are presented in 
Table [8] Obtaining the same accuracy as for the ground state sometimes requires 
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Table 7: Pressure acting on a helium atom confined in a sphere of radius R in a.u. and 
atm. The powers of ten are indicated between brackets. 


R 

N 


P (a.u.) 

P (atm) 

Ref. [in] 

Ref. [H 

Ref. [m 

0.1 

20 

20 

1.507426 73864(5] 

4.377006 84319(13] 




0.2 

20 

20 

4.512 880 064 60(3] 

1.310 372 59848(12] 




0.3 

20 

20 

5.682 940 282 5 (2] 

1.6501145872 (11] 




0.4 

20 

20 

1.287118 2851 [2] 

3.737312 996 8 [10] 




0.5 

20 

20 

4.017215 609 5 [1] 

1.166450 0677 [10] 

1.16640(10] 

1.167(10] 


0.6 

20 

20 

1.534428 7331 [1] 

4.455410 5969 [9] 

4.45436(9] 

4.455[9] 


0.7 

20 

20 

6.732 316 374 6 

1.954 814392 4 [9] 

1.954 50(9] 

1.954[9] 


0.8 

20 

20 

3.266 880 6106 

9.485 806 7871 [8] 

9.48410(8] 

9.485[8] 


0.9 

20 

20 

1.7109177385 

4.967869 056 3 [8] 

4.96703(8] 

4.967[8] 


1.0 

20 

20 

9.510 085 6621 (-1] 

2.761375 326 3 [8] 

2.76049(8] 

2.762 [8] 


2.0 

25 

25 

1.383499 912 (-2] 

4.017169411 [6] 

4.014 56(6] 

4.018[6] 

4.00729(6] 

3.0 

25 

25 

6.192 83139 (-4] 

1.798168 01 [5] 

1.79150(5] 

1.798 [5] 

1.793 91(5] 

4.0 

25 

25 

3.702 24159 [-5] 

1.074 993 32 [4] 

1.086 77(4] 

1.074[4] 


5.0 

30 

30 

2.365111 17 [-6] 

6.867403 58 [2] 

9.02194(2] 

6.882[2] 

6.850 3 [2] 

6.0 

30 

30 

1.526 228 5 (-7] 

4.4316001 [1] 

1.33015[2] 

4.434[2] 


7.0 

35 

35 

9.871147 (-9] 

2.866 214 



2.84795 

8.0 

35 

35 

6.399 98 (-10] 

1.858 31 [-1] 


1.868(-1] 


9.0 

35 

35 

4.1614 (-11] 

1.2083 [-2] 




10.0 

35 

35 

2.711 (-12] 

7.874 [-4] 





higher numbers of mesh points. At R = 0.1, the excitation energy exceeds 1000 a.u. 
For R = 3, 5, 7 and 10, a comparison is possible with Ref. [16]. The six-digits results of 
Ref. [T6| agree with the present ones at the level of 10“®. The best variational energies 
of Ref. [I3| have an accuracy around one percent. 

Energies and mean distances for the lowest triplet level are gathered in Table O 
They are obtained from the lowest-energy solution for the spatially antisymmetric state. 
The numbers of mesh points are very close to those for the ground state. At R = 0.1, 
the excitation energy reaches about 1500 a.u. The 2^S state is thus above the 2^5 
contrary to the free-atom situation as already observed in Ref. [T3|. This is probably 
due to the additional constraint that the wave function must vanish for ri = r 2 in a tiny 
space. The usual ordering is recovered around R = 1 with an excitation energy of about 
5 a.u. At all confinements, the triplet state is slightly less extended than the singlet 
excited state like for the free atom [27| . The best variational energies of Ref. [13] reach 
an accuracy of abour 10“^ above R = 1. Below that value the perturbation approach 
is more accurate; its relative accuracy improves from about 10“^ at R = 1 to about 
10-^ at R = 0.1. 

7 Concluding remarks 

Various types of helium confinement can be accurately treated with the Lagrange-mesh 
method. They are based on the perimetric coordinates for soft confinements or on a 
new coordinate system for the hard confinement. The method rapidly leads to a large 


16 





Table 8: Energy and mean interparticle distances of the 2^5* singlet level of a helium 
atom conhned in a sphere of radius R. 


R 

N 


E 

(G2) 

(ri) = (rs) 

0.1 

15 

15 

1963.757922 04146 

0.077448 989 286 0 

0.059 004994435 03 

0.2 

15 

20 

477.023 366 387052 

0.155 0472684474 

0.117679 576 745 69 

0.3 

15 

20 

205.784321398 280 

0.232 963 5031645 

0.17601138786641 

0.4 

15 

20 

112.226 705 249 011 

0.311436 535 890 

0.233 983451406 40 

0.5 

15 

20 

69.550421013 503 

0.390 828 529 869 

0.291568 743 866 7 

0.6 

15 

20 

46.705 530 340 508 9 

0.471731678 057 

0.348 716 226 254 7 

0.7 

15 

20 

33.131390 885 993 4 

0.555 200 036 019 

0.405 310 225 626 3 

0.8 

15 

20 

24.447 862 429155 2 

0.643 293 324656 

0.461031786 9316 

0.9 

15 

20 

18.574 582 880 828 5 

0.740 350 389 720 

0.514 775 812 4513 

1.0 

15 

20 

14.413 766 091552 3 

0.853 808 811249 

0.561631072 086 

2.0 

20 

25 

0.946 588473 886 9 

1.343 903122 341 

0.883 824224122 

3.0 

25 

25 

-1.114121513 7421 

1.866 351750 214 

1.203 075 085 099 

4.0 

25 

25 

-1.717517590 2269 

2.408 986 509 893 

1.498 595 969 032 

5.0 

30 

30 

-1.949 7613214615 

2.948 630420 83 

1.782 952 066 008 

6.0 

30 

30 

-2.050 702179185 3 

3.435 878 96714 

2.036478158 792 

7.0 

30 

30 

-2.098 085 0744321 

3.858190 026 36 

2.254322 411040 

8.0 

30 

30 

-2.121511222 458 8 

4.215 842 573 24 

2.437665 074483 

9.0 

35 

35 

-2.133 453 292 480 3 

4.510 57216717 

2.588 081892182 

10.0 

35 

35 

-2.139 619 886178 6 

4.744 712 44944 

2.707196 655 205 


but sparse symmetric matrix. The eigenvalue problem is not generalized, even with the 
new coordinate system for which the basis is not orthogonal. This striking property 
results from the systematic use of the Gauss quadrature associated with the Lagrange 
basis. This simplifying approximation does not cost accuracy mu- 

The present results improve previous works by at least hve orders of magnitude for 
the energies and allow a very accurate calculation of mean interparticle distances. This 
is obtained with short computer times on a personal computer. The scale parameters 
given in the present tables allow avoiding a search for their optimal domain of values 
in future similar calculations. Together with the energies, approximate wave functions 
are also available. Their high accuracy is conhrmed by the high accuracy on the mean 
values of the coordinates. They are available for other applications. 

A serious challenge is the extension of the method to more than three charged 
particles because a convenient coordinate system, i.e. where the different Coulomb 
singularities can be regularized, is not available yet. 
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Table 9: Energy and mean interparticle distances of the 2^S triplet level of a helium 
atom conhned in a sphere of radius R. 


R 

N 


E 

(C2) 

(ri) = (rs) 

0.1 

15 

15 

2370.727023 303 97 

0.074235 935 806 89 

0.049 746 873483 20 

0.2 

15 

15 

568.187610 564 046 

0.148106 984196 29 

0.098 978 62549715 

0.3 

15 

15 

241.491539 024 040 

0.221580 296 180 99 

0.147681956 05843 

0.4 

15 

15 

129.543 621 160 698 

0.294623 610 98139 

0.195 843 625489 01 

0.5 

15 

15 

78.820 863 703 3835 

0.367205 348 9194 

0.243 450 531053 3 

0.6 

15 

15 

51.856 719 933 0814 

0.439 294 7116545 

0.290 489 790185 5 

0.7 

15 

15 

35.951107383 7302 

0.510 861789 0944 

0.336 948 829 8349 

0.8 

15 

15 

25.855 504 0261363 

0.581877672 009 0 

0.382 815 481229 5 

0.9 

15 

15 

19.089 233 032 7807 

0.652 314 569 086 5 

0.428 078 0791261 

1.0 

15 

15 

14.359 714 920 8699 

0.722145 926 846 8 

0.472 725 564359 7 

2.0 

15 

20 

0.560 251233 73470 

1.382427667612 2 

0.883 572 145 6814 

3.0 

20 

25 

-1.370 510 61013706 

1.965 771 607822 

1.228 588 674995 6 

4.0 

20 

30 

-1.874611596 37791 

2.474178173 270 

1.515 265 802 863 

5.0 

20 

30 

-2.048 044 093 51074 

2.916 562 529 64 

1.755 588 974154 

6.0 

30 

30 

-2.117816 28516207 

3.298 089 015 33 

1.957801932 623 

7.0 

30 

30 

-2.148 564470 05490 

3.619 012 948 88 

2.125 320 088136 

8.0 

30 

30 

-2.162 783 826 65888 

3.878 261844 64 

2.259 362 209113 

9.0 

30 

30 

-2.169481052 37618 

4.076 816120 79 

2.361388 47941 

10.0 

30 

30 

-2.172 627252 24936 

4.219 554 50184 

2.434424170 24 
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